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Abstract. We present H and K band imaging spec- 
troscopy of the core regions of the cD/AGN galaxy 
NGC 1275. The spectra, including hues from Hs, H, ^^CO 
bandheads, [Fell], and [Feiii], are exploited to constrain 
the star formation and excitation mechanisms in the 
galaxy's nucleus. The near-infrared properties can largely 
be accounted for by ionized gas in the NLR, dense molec- 
ular gas, and hot dust concentrated near the active nu- 
cleus of NGC 1275. The strong and compact H2 emission 
is mostly from circumnuclear gas excited by the AGN and 
not from the cooling flow. The extended emission of late- 
type stars is diluted in the center by the thermal emission 
of hot dust. The line ratios [Fe Ii]/Br7, as well as [Fe ii]/H2 
can be explained with X-ray excitation from the central 
AGN. 

Key words: Galaxies: Individual NGC 1275 — Galaxies: 
Cooling Flows — Galaxies: Active — Infrared: Galaxies 
— Infrared: ISM: Lines and Bands — Instrumentation: 
Spectrographs 



1. Introduction 

NGC 1275 (Perseus A, 3C 84) at the center of the Perseus 
cluster has been variously classified as an active galaxy, 
a merger, a cooling flow galaxy, a BL Lac object, and 
an FRI radio jet source. However, many of its phys- 
ical characteristics — from X-ray to radio — remain 
baffling. NGC 1275 e xhibit s the characteristics of several 
source type s: Sey fert (|1943|) included it in his original list, 
Minkowski ( 1957 ) called it a galaxy "collision". The dom- 
inant low velocity (LV) component (w 5200 km s~^) is the 
giant cD elliptical NGC 1275, the other component (HV, 
« 8300 km s^^) is probably a nearby spiral galaxy (at 
lea st par tially) in front of the main galaxy (Haschick et 
al. |1982D . Veron ( |1978[ ) suggested that NGC 1275 was a 
BLLac object, while other authors have suggested that it 
was more compatible with LINERs (e.g. Norgaard-Nielsen 
et al. |1990[). 



NGC 1275 has a powerful, nearly point like X-ray emit- 
ting nucleus with a flux of Lj^ « 2x10"*^ erg s^^ (Prieto 



1996), which is embedded in an extended halo of diffuse 



X-ray emitting intracluster gas. The nucleus of NGC 1275 
is displaced by about 25" from the center of symmetry 
of the large scale 2') cluster gas halo, perhaps due 
to its peculiar motion through the gas (Bohringer et al. 
1993). Temperature profiles of the X-ray halo suggest that 
a "cooling flow" of « 200MQyr~^ is condensing out 
of the surrounding intracluster medium onto the galaxy 
(Fabian et al. |1981|, Cardiel et al. |1995[). As in aU cool- 



1, Cardiel et al. 
ing flow galaxies, the ultimate fate of this purported in- 
flowing material is unclear, but visible band colors and 
Balmer absorption lines spectra characteristic of AO stars 
show that at most 10% of the gas ca n for m stars with a 
standard Salpeter IMF (Romanishin |l986| ). NGC 1275 is 
unique among cooling flow galaxies in that it in fact shows 
evidence of ongoing star formation. 

Even by the standards of infrared luminous mergers, 
NGC 1275 has very powerful, extended II2 emission (Fis- 
cher e t al. |1987| ) matched only by NGC 6240 (van der Werf 
et al. 1993| ). In NGC 1275, this emission has been specu- 
lated to have its origin in the cooling flow, although fast 
J-shocks, UV pumping, and X-ray heating via supernovae 
and/or AGN also produce strong II2 emission. Merging 
events are a likely mechanism for supplying NGC 1275 
with substantial (M >, 3x10^Mq) quantities of II2 in 
giant molecular clouds which are detected through mm 
wave CO observations within the central regions (Inoue 
et al. 1996). CO data also show two kinds of motion: a 
large ( <L 5 kpc) scale rotation which has been attributed 
to merging/cannibalism, and smaller scale turbulence, at- 
tributed to the inflow of cooling flow gas (Renter et al. 
|1993D . Recent 2.6 mm CO interferometric observations re- 
veal the majority of the gas to be in a ring structure of 
radius « 4" oriented in the East- West direction. The ring 
has a westerly extension up to 30" from the nucleus, and 
nearly 80% of the total CO is on the western side (Inoue 



et al. 1996). Very deep interferometric observations of the 
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nucleus itself show no CO absorption from very cold gas, 
which puts very tight limits on any mass deposition from 
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Table 1. Observations of NGC 1275 with the MPE 3D 
imaging spectrograph. 



fiir Astronomic at Calar Alto, Spain. The CHARM in- 



Date 


Band 


*int 
sec 


Seein 


15/01/95 


K 


2400 


1'.'4 


15/01/95 


H 


1100 


1'.'5 


16/01/95 


K 


2000 


1'.'5 


16/01/95 


H 


1200 


1'.'6 


21/01/95 


H 


2400 


I'.'l 



Calibrator 



my 



Mel 20497 (F6V) 6.9 

Mel 20497 (F6V) 6.9 

PPM 45428 (F8V) 4.1 

Mel 20497 (F6V) 6.9 

Mel 20497 (F6V) 6.9 



a cooling flow if the gas cools to nearly 3 K (Braine et al. 



1995) 



NGC 1275 is a Fanaroff-Riley type I (F-RI) radio 
source with powerful (« 0.4 Jy at 22 cm) radio lobes 
lOkpc long and oriented at position angle (P A) 160 ° 



19901) . 



which emanate from a central core (Pedlar et al. 
The asymmetry of the radio lobes and their spectral in- 
dices suggest that we are looking nearly pole- on at the 
centra l F-R I radio engine (Pedlar et al. 199C , Levinson 
et al. 1995 ). Within the central 30" region, relativistic 
particles and the magnetic field from the bipolar radio 
lobes produce a pressure, which exceeds the thermal gas 
pressure and so must disturb any cooling flow within that 
region (Pedlar et al. 1990 , Bohringer et al. 1993 ). At even 
smaller scales, VLBI imaging shows that the central core 
(w 1") breaks up into many smaller components and a 
southern/northern jet pair (Venturi et al. 1993, Vermeulen 
et al. 1994). Recently, a new component has been dis- 
covered only 55 mas south of the central peak (Taylor & 
Vermeulen 1996). Near infrared (NIR) annular photome- 
try shows that the nuclear K-band light is dominated by 
non stellar emission which could be either thermal emis- 
sion from hot dust, or optically thin synchrotron radiation 



(Longmore et al. 1984). Recent HST observations in the 
UV show that any nuclear engine/star cluster must be 
<17pc in size (Maoz et al. |1995D . 

Our goal in this work is to better constrain the many 
existing models of NGC 1275 with near infrared (NIR) 
imaging spectroscopy and spectral synthesis modeling. In 
dusty environments, NIR imaging has optical depth pene- 
tration 10 times deeper than visible wavelengths, and NIR 
spectroscopy is sensitive to a wealth of diagnostic lines 
which describe the physical state of the gas and stars in 
the nucleus. Throughout this paper we assume a scale of 
344pc"-i (Ho = 75kms-iMpc^i, z = 0.0172). 

2. Observations and data reduction 

We observed NGC1275 using the Max-Planck-Institut fiir 
extraterrestrische Physik (MPE) imaging spectrometer 
3D and the Max-Planck-Institut fiir Astronomie (MPIA) 
tip-tilt image tracker "CHARM" on 1995 January 15, 16 
and 21 at the 3.5m telescope of the Max-Planck-Institut 



strument (McCaughrean et al. 1994) removes atmospher- 
ically induced image motion and produces stable images 
of increased spatial resolution. We locked CHARM on the 
bright point-like nucleus of NGC1275, and corrected at 
5 Hz. Under seeing conditions varying from 1.0 to 1'.'7, we 
Achieve a final overall image quality of FWHM 1'.'5 as de- 
termined by fitting to the point-like nucleus (Maoz et al. 
1995|). The "3D" instrument (Krabbe et al. |1995| , Weitzel 



et al. 199£) is an field imaging spectrograph which pro- 
vides, during in a single exposure, H or K band data cubes 
with spectral resolution of R = A/AA = 1000 and 2000. 
The two spatial cube dimensions are 16 pixels on a side, 
with pixel sizes selected to be 0.5" (giving an 8" x 8" field 
of view), while the third dimension of 256 pixels is spec- 
tral. The intrinsic spectral sampling of the camera is with 
pixels of size A/R, so by dithering the spectral sampling by 
1/2 pixel on alternate data sets, we achieve fully Nyquist 
sampled spectra. The spectral resolution, achieved with 
the current reduction procedure, is 670, or 450km s~^; 
future software will produce full R 1000 spectra after 
co-addition. For extended objects, the efficiency gain over 
a long slit spectrograph is substantial, as at any point 
within our 8"x8" field, it is possible to obtain the full 
band object spectrum. 

The data (Table ^ frames were interleaved with iden- 
tical sky exposures 60" away to the West. The source 
observations were preceded and/or followed immediately 
by a set of spectral calibration observations of nearly fea- 
tureless stars of known brightness. The few known spectral 
features in the reference star spectra (e.g. Br7 in K band 
and Br9 - Brll in H band) were later removed by interpo- 
lation. A temporally variable deep CO2 absorption near 
2.01 /im made accurate sky subtraction in this regime dif- 
ficult. In addition, the temporal variability of OH lines in 
the H band produces spurious sky noise which dominates 
the background across the band. 

A run-down on typical 3D data reduction with the 
Groningen GIPSY software is given in Weitzel et al. 
( 1996| ). Spectral baselines were calibrated using the color 
temperature of the calibrator stars and a featureless 
Nernst glower. The absolute wavelength calibration is ac- 
curate to better than 1/4 pixel. Individual data sets were 
properly weighted, registered onto the central peak of 
NGC 1275, and co-added. 

The K and H band spectra of the nucleus within a 3" 
aperture are presented in Figs. |^ and|^, which present the 
enormous variety of emission and absorption lines. The 
signal to noise in the K-band spectrum is « 70 in the con- 
tinuum. For the H band, the signal to noise is ~ 25 in the 
continuum. 



3. Nature of the continuum emission 

Fig. ^ represents the true, line emission free, K band con- 
tinuum structure of NGC 1275. Its peaking at the nucleus 
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Fig. 1. K band spectrum of the cen- 
tral 3" of NGC 1275. This aperture 
contains nearly all of the line flux. 
A linear baseline continuum of mean 
value 9.0 x 10"^'' ergs~^ cm"^ A"^ 
has been removed. 




Fig. 2. H band spectrum of the cen- 
tral 3" of NGC 1275. This aperture 
contains nearly all of the line flux. 
A linear baseline continuum of mean 
value 1.3 x 10~^^ erg s~^ cm^^ 
has been removed. 

by a power-law emissivity of the form e(j^g^ oc A^^. From 
these contributions a synthetic galaxy spectrum including 
extinction effects has been modelled and (least square) 
fitted to the observed continuum spectrum. Details of the 
procedure are given in the appendix, contrains on the pa- 
rameters are bein g d iscussed below, and the results are 
presented in Sect.|3.2|. 

3.1. Physical constrains on the continuum emission 

The size of the search space was reduced by placing con- 
strains on the total extinction, Galactic extinction, free- 
free contribution, dust temperature, and stellar popula- 
tions. 

i) The Galactic extinction toward NGC 1275 is Ag 
= 0.7, or Ay = 0.53, and there is an additional com- 
ponent of extinction intrinsic to NGC 1275 itself. From 



confirms that the bulk of the emission arises in the cen- 
tral 300 pc. The slight south-east north-west extension. 



almost along the radio PA of 160° (Pedlar et al. |1990D , is 
2(7 significant. 

We analyzed the the continuum spectrum with respect 
to contribution from stars, the active galactic nucleus, hot 
dust from the circumnuclear environment, and free- free 
emission. K and M stars (Tg-g? k, 4000) , dominating the K 
band stellar emission, show strong CO absorption features 
in their more or less BB spectrum which can be used to 
constrain the fraction of total light originating in stars. 
NIR spectra from AGN show a power law spectrum of the 
form 1\ oc A" where a varies between —0.5 and —1.5 (Laor 



& Draine |1993| , Acosta-Pulido et al. |l99g). Hot dust in this 
context is mainly referred to as being reprocessed radia- 
tion from AGN and, to a minor extent, from starburst 
activity. It's spectrum is that of a black body multiplied 
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Table 2. Line strengths in the central 3" aperture 



Species 


A 


a 


Width 


h 




Fell 1.534 


1 


5335 


660 


1.3 


0.5 


Fell 1.644 


1 


6439 


526*^ 


63 


2.5 


Ha 1-0 S(7) 


1 


7480 


256 


11 


0.7 


BrS 


1 


9446 


263 


3.1 


0.5 


Ha 1-0 S(3) 


1 


9576 


301 


42.5 


2.5 


Ha 1-0 S(2) 


2 


0338 


225 


14.4 


0.7 


Hel 


2 


0584 


519 


2.8 


0.5 


Ha 2-1 S(3) 


2 


0735 


183 


3.1 


0.6 


Ha 1-0 S(l) 


2 


1219 


238 


41.1 


2.5 


Ha 2-1 S(2) 


2 


1540 


517 


3.4 


0.5 


Br7 


2 


1661 


982 


9.7 


2.5 


Ha 3-2 S(3) 


2 


008 


207 


1.3 


0.5 


Ha 1-0 S(0) 


2 


2231 


328 


10.0 


0.7 


Felll 2.244 


2 


243 


504 


2.2 


0.5 


Ha 2-1 S(l) 


2 


2473 


754 


7.7 


0.6 



^ fim observed corrected for z — 0.0172 

^ km s~^ corrected for instrument resolution of 450 km s~^. 
Width error typically 30 - 50 km s~^ 



10" 



erg s cm 



Best fit of single Gaussian. Two component fit of 1345 ± 



120 km s and 280 ± 22 km s is superior. See text. 



UV measurements we have E(B-V) 

1995| ) 



1995, Levinson et al. 



0.2 - 0.3 (Maoz et 
which both suggest that 



Ay « 1.0 at the nucleus. Integrated over the central stel- 
lar area, H -K is 0.15 redder than typical ellipticals (Long- 
more et al. 1984 ), which requires Ay w 2.4. Hence we con- 



strain the V band extin ction to Ay < 2.0, or Aj^ < 0.2 



(Rieke & Lcbofsky 1985 ). Since we know that the central 
lig ht is l ikely to be dominated by the AGN (Longmore et 
al. |1984D , we choose to use a screen model of extinction 
there, while farther out at distances 500 pc we will use 
a mixed stars/gas model. 

ii) The free-free continuum can be obtained from the 
Br7 line flux using case B recombination and assuming 
T = 10000 K. (Joy & Lester |1988D have shown that this 
method gives consistent results compared to a NIR multi- 
color continuum decomposition (Joy & Harvey 1987) and 
radio measurements (Wynn- Williams & Becklin 1986| ). 
Using the Br7 line flux from Table |^ we can thus constrain 
the free-free contribution to be about 2% of the underly- 
ing continuum. There is some suggestion that the free-free 
emission could be higher in NGC 1275 because the central 
plasma is constantly being replenished with free electrons 
via the cooling flow gas. However, if there were a substan- 
tial population of such ions, they would contribute to the 
scattering of the central source at the 1% level (Sarazin & 
Wise 1993), which would in turn produce a blue excess. If 
the 1% scattering is responsible for the observed nuclear 
blue excess of ~ 10^^ ergs~^ (Crawford & Fabian 1993D , 



then the central source would require an X-ray/UV lumi- 
nosity of roughly 10"'^ergs~^, which is « 600 times the 




Fig. 3. The K band continuum of NGC 1275. Orientation: 
N up, E left. FWHM 1'.'6. Field of view 8"x8". Contours 
are at 1.25, 2.5, 5, 10, 20, 30. . . % of peak flux of 8.5 x lO^i'"^ 
erg s~^ cm~^ arcsec"^. 



point X-ray luminosity. Hence we set the free-free emis- 
sion to a typical value of < 3% of the total continuum, or 
sometimes ignore it completely, as our fits constrain the 
relative contributions to only w ±10%. Had we not done 
so, the free-free and AGN emission would have traded off 
influence with each other, as they are both power laws of 
similar slopes. 

in) Hot dust has turned out to be very common in the 
nuclear regions of AGNs (Granato et al. 1997, Maiolino 
et al. 1995) and, to some extent, in starburst galaxies as 
well (see e.g. Boker et al. 1998] ). It surrounds the nucleus 
forming a torus which absorbs radiation from the AGN 
and reemits it in the mid-infrared. Since hot dust too close 
to the AGN will start to sublime, the maximum dust tem- 
perature is set by the sublimation temperature of the dust 
material. Silicate and graphite grains, which ar e the main 
constituents of interstellar dust (Laor & Draine 1993 ), will 
sublime at temperatures of about 1400K and 1750K resp., 
setting the limit for the the dust temperature at about 
1500K at the inner edge of the dust torus. The minimum 

1 /2 

radius for dust grains to survive is ^ 0.2L4g pc (Laor & 
Draine 1993D , where L46 is the central source bolometric 



luminosity in units of lO''^ ergs ^. For NGC1275 with 



J46 



0.04 (Levinson et al. 1995), this is > 0.05 pc 



The dust temperature observed will be a mixture of dust 
at different temperatures along the line of sight weighted 
with their column density distribution. Thus the effective 
temperature will be lower and, e.g. for NGC 1068, in the 
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range of 700 K (Thatte et al. |1997| ). As NGC 1275 is a 
F-RI type radio source, we are probably looking nearly 
pole on at the central source, and hence see the largest 
possible fraction of the hot inner surfaces of the clouds 
which make up an obscuring torus (Pier & Krohk 1992D , 
suggesting that we limit the temperature of the dust in 
our simulations to a relatively high value of 1000 K. 

iv) The NIR flux of essentially all integrated galaxy 
populations is dominated by gi ants a nd, to a lesser ex- 
tent, by supergiants (Oliva et al. 1995 ). We rule out K5 V 



stars because for all reasonable stellar populations they 
are too faint to contribute substantially to the overall H 
and K band fluxes. As an underlying stellar population we 
choose K5 III stars because they dominate, to first order, 
the light of ne arly all galaxies at H and K bands (Fro- 
gel et al. |1978|, Oliva et al. 



1995) 



This is due to the fact 
that their effective temperatures of w 4000 K have a black 
body emission peak in those bands, and because they are 
typically 200 times brighter than K5 V stars with similar 
effective temperatures. Digitized libraries of NIR stellar 

spectra are available in the K band 1986, so we used the 

STAT? 

spectrum of the K5 III star 7 Dra as our template IJ 

in the K band. 

In the H band, no such digitized stellar library was 
publicly available at the time of this data reduction, so 
we constructed a continuum spectrum appropriate to the 
stellar population of NGC 1275 by summing the contribu- 
tions of many black body curves of different temperature 
scaled by the fractional number of stars having that tem- 
perature. This produces a continuum spectrum without 
any lines. The H band CO absorption features are rela- 
tively weak, so they would have contributed little to the 
fitting anyway. The galaxy population, necessary for this 
procedure, was modeled with the stellar population evo- 
lution code STARS (Sternberg & Kovo, in preparation). 
Given a set of input parameters it determines the binned 
number density function of stars in the HR-diagram. Sum- 
ming the series of appropriately scaled black body curves 
at the various temperatures results in a H-band contin- 
uum that approximates the underlying H band continuum 
of NGC 1275, if the input parameters to STARS are varied 
adequately. The composite spectrum is then normalized to 
the same wavelength as the K band stellar library spec- 
trum in order to produce the final H and K band underly- 
ing stellar spectrum, IJ -l^j^ All K spectra are normalized 
to a wavelength of 2.27/im, and all H spectra are normal- 
ized to 1.57/xm, which we chose because the spectra are 
free of lines there. 



3.2. Decomposing the continuum emission 

In the central 1" (0.3 kpc) region we expect that the con- 
tribution of the AGN overwhelms the the stellar popula- 
tion. This is indeed the the continuum emission 
from the center of NGC 1275 is well modeled by the light of 




Fig. 4. Spectral fitting to the NGC 1275 spectrum in an 
aperture of 1" centered on the nucleus. Best fit is achieved 
with a relative stellar contribution of 0.0, an AGN power- 
law (see text) contribution of 0.5, a hot dust contribution 
of 0.5 at 700 K, and a K-band extinction of 0.05 (screen 
model) . The plot shows (in decreasing order from the top) 
a) the raw spectral data; b) the final derived model fit to 
the continuum; c) the statistical weights given to each 
data point; d) the derived power-law component of the 
model; e) the derived dust component of the model; f) 
the residuals between model and data. Note that because 
we are fitting the continuum, not the lines, the weights in 
regions of strong line emission have been set to zero (cf 
Eq. (A4)). The spectra are all normalized at A = 2.27/im, 
so the relative contributions of the various components 
can be read off at that wavelength. 



K5 III or K5 lb stars heavily diluted by emission from hot 
dust and AGN typical power-law emission (Fig. |4|). The 
best fitting values are cg^p^j^ = O.OiO.l, cp^ = 0.5±0.1, 
CJJJ3 = 0.5 ± 0.1, = 700 ± 200, = 0.05 ± 0.1, where 
the Cj represent fractional contributions and were the er- 
rors are those derived from the estimated systematic er- 
rors in the spectra and from multiple trials with slightly 
varied weights and initial guess parameters. The relative 
insignificance of spectral features in the underlying stellar 
population shows that the fitting process at the nucleus 
is dominated by the power laws, and that we are unable 
to distinguish between the various possibilities of under- 
lying populations. The relative contribution of power-law 
emission is not particularly sensitive to the adopted AGN 
spectral index (-0.5 through -1.5, see appendix), because 
its main function is to dilute the very steep slope of the 
stellar continuum (very nearly a black body). 

In contrast, the same fitting process applied to a 1" 
wide annulus 3" (0.9 kpc) from the nucleus yields com- 
pletely different result. Because the CO absorption band- 
heads are clearly observed in the K band spectrum (they 
are much less diluted), we fit this region only using the K 
band data. The fit is excellent (Fig. g) , and yields the fol- 
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Fig. 5. Spectral fitting to the NGC 1275 speetrum in a 1" 
wide annulus of radius 3" centered on the nucleus. Best 
fit is achieved with a relative stellar contribution of 0.6, 
an AGN power-law (see text) contribution of 0.4, a hot 
dust contribution of 0.0, and a K-band extinction of 0.01 
(mixed stars and gas model). 



lowing contributions: cgrp^j^ = 0.6±0.1, cp^ = 0.4±0.1, 
cjjQ = 0.0 ± 0.1, Tjj ~ 0.01 ± 0.1, where we have cho- 
sen to apply a mixed stars/gas model of extinction rather 
than a simple screen. Again, we assume a K5 III star 
dominated underlying stellar population. Note that the 
AGN contribution could also be interpreted as free-free 
emission at this distance from the nucleus, where we do 
not expect any direct or scattered contribution from the 
AGN. By the time we reach 4" (1.2 kpc) from the nucleus, 
the stellar component dominates completely, and we find 
CSTAR = 0-90 ± 0-l> = 0.1 ± 0.1, cjjp, = 0.0 ± 0.1, 
rj^ = 0.01 ±0.1. Hence the behavior of the the various 
components is as expected from a central AGN dominated 
source. 

4. The H2 emission 

The H2 emiss ion in NGC 1275 is exceptionally luminous 
(Fischer et al. 1987) even if compared with other powerful 
H2 emitting galaxies. The excitation of H2 in AGNs is 
usually explained to be a mixture of shocks (Brand et al. 
1989| ), UV photons (S ternb erg & Dalgarno |l989| ), and X- 



rays (Draine & Woods |1991| ) . The shocks can be caused by 
galaxy merging or by supernovae, the UV radiation could 
come from OB star associations (Puxley et al. 1990), while 
only an AGN with a black hole can produce the quantities 
of X-ray radiation needed to excite the H2 seen in Seyfert 
nuclei. 

J^.l. Spatial distribution of H2 emission 

The H2 line emission (Fig. ^ is strongly dominated by 
an unresolved source centered on the nucleus surrounded 
by weak extended emission extending a few arcseconds 



Fig. 6. Line map of the H2 1-0 S(l) transition. The emis- 
sion is almost completely symmetric, with FWHM of 1'.'6. 
The other H2 transitions show similar structure and ex- 
tent. Contours are at 2.5, 5, 10, 20, 30. . . % of peak flux 
of 2.1 X 10^^'* erg s~^ cm~^ arcsec"^. 



west. The fact that the primary exciting source is not dis- 
tributed over scales 300 pc suggests that most of the 
H2 emission is associated with the circumnuclear environ- 
ment and probably the active galactic nucleus and not, as 
speculated before (Inoue et al. 1996 ) with the coohng flow. 
The low level extended emission reach ss 2% of the peak 
flux at our limiting distance of 5" or 3.5 kpc (Fig. The 
total H and K band H2 emission amounts to 1.3 x 10~^^ 
ergs~^ cm~^ (see Table p|) or, for a distance of 70Mpc, a 
total luminosity of Ljj_^ = 



'S X 10"*" ergs-i. 



4-2. Excitation of the H2 emission 

The high signal-to-noise in our spectra allow us to check 
the effective excitation temperature from the column den- 
sity of molecules in each quantum state (see e.g. Beckwith 
et al. Il97|): 



Nr.. 



Anl/Auihu 



(1) 



where / is the surface brightness in the line and A^j^ is 
the Einstein A coefficient. The surface brightnesses have 
been dereddened using Ay = 2.0 (Sect. ^) as a realistic 
value. Screen and mixed extinction models (McLeod et 



al. 1993) for the dereddening gave flux correction between 



10% and 20%. Since the appropriate model was uncertain, 
the correction was averaged between screen and mixed 
model and the uncertainties propagated in the errors. 
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Table 3. Column density and temperature analysis of H2 



Line 


Aa 




'"(^col/ffu) 


Tex 


1-0 S(0) 


2.5(-7) 


3.8(15) 


43.47 


6472 


1-0 S(l) 


3.5(-7) 


1.1(16) 


43.09 


6951 


1-0 S(2) 


4.0(-7) 


3.2(15) 


42.72 


7585 


1-0 S(3) 


4.2(-7) 


8.7(15) 


42.42 


8365 


1-0 S(7) 


3.0(-7) 


2.9(15) 


40.79 


12817 


2-1 S(l) 


5.0(-7) 


1.5(15) 


41.10 


12551 


2-1 S(2) 


5.6(-7) 


5.6(14) 


40.98 


13150 


2-1 S(3) 


5.8(-7) 


4.9(14) 


39.53 


13887 


3-2 S(3) 


5.6(-7) 


2.0(14) 


38.66 


19089 



^ Turner et al. 



1977 



H2 flux dereddened with Ay = 2.0 



Fig. 7. Azimuthally averaged profile of the H2 1-0 S(l) 
emission. There is a very low level extended component 
to the emission, which reaches a few percent of the peak 
flux at r = 4". 



In thermal equilibrium, the line ratios should be given 
solely by the number of atoms in each quantum state and 
their spontaneous emission coefficients if the H2 is op- 
tically thin. From Table we conclude that this is in- 
deed the case. The column densities are always lower than 
N^,oi(H2 1-0 S(l))= 9.6 X 10^5 cm-2. Summing over all the 

If 



observed II2 lines gives N(,qj(II2^°^) « 3 x 10^^ cm 
we assume that this gas is uniformly distributed within 
the emitting region, then the total mass of contributing 
H2 would be - 1000 Mq. 

In Fig. H we plot the natural log of N^^^j per mode 
against the energy of the upper level giv en as a tempera- 
ture (compare with Kawara & Taniguchi 1993). The error 
bars include the errors given in Table |l and the dered- 
dening. The inverse slope is the effective excitation tem- 
perature Tex of the gas and corresponds to the kinetic 
temperature if local thermal equilibrium (LTE) prevails. 

Depending on the conditions within a molecular cloud, 
the initial ortho/para ratio of 3 of molecular hydrogen at 
formation may evolve and eventually drop to lower val- 



ues (Martini et al. 1997). Exchange reactions with H and 
H+ (Flower & Watts |1984| ) or a low temperature environ- 
ment (T ^ 200 K; Burton et al. 1992) can even lower the 
ratio to below 2.0 (see e.g. Hora & Latter |1996| ). A high 
temperature environment in turn will keep the ortho/para 
ratio at 3. In a photodissociation region (PDR), however, 
freshly exposed H2 m ay kee p the ortho/para ration down 
(Chrysostomou et al. 1993| ) although the environment is 
hot. In Fig. 1^ the column densities have been divided by 



the statistical weights corresponding to their states which 
include the spin degeneracies gs = 1 for even (para) and 
gs = 3 for odd (ortho) J. If the ortho/para ratio in the 
molecular hydrogen is 3:1 then the ortho and para column 
densities should lie along a line. For the 1-0 transitions, 
were we have enough data points, this is indeed the case. 
This demonstrates that the either the molecular hydrogen 
has just been formed or that the II2 has been kept in a 
high temperature environment. 

The distribution of data in Fig. || shows that the H2 
gas is nearly isothermal. Below 8000 K upper level energy, 
the excitation temperature is Tex ~ 1500 K, but above 
8000 K, Tex does only vary between 2600 K and « 2900 
K. If no extinction correction is applied, the temperatures 
come out to be 1450 K, 2400 K, and 2650 K resp., which is 
still in agreement with the la errors given in Fig. ||. All 
date (except for 2-1 S (3), see below) out to level energies 
of « 20000 K lie along an only slightly curved line. In par- 
ticular, the column density per mode of 1-0 S(7) is close 
to those of 2-1 S(l) and 2-1 S(2), indicating that the ro- 
tational and vibrational temperatures are about the same 
and the fluorescent excitation of the gas is rather small 
(Draine & Bertoldi 1996). The proposed likely mechanism 
to populate the levels of this hot II2 gas with Tex < 2700 
K is coUisional excitation in a post shock environment, 
which has been observed in many galactic sources (see 
e.g. Fischer et al. |1987|, Gredel et al. |1994 Wright et al. 



1996). However, at column densities prevailing in the cen- 
tral 300 pc of NGC 1275, a PDR scenario cannot be com- 
pletely excluded. In dense PDRs, the vibrational ladder 
is thermalized while the rotational ladder in a given vi- 
brational state may be out of thermal equilibrium (see 



e.g. Timmermann et al. 1996). The fact that we do not 



see gas of temperature above 5000 K (Wright et al. 1996) 



may, however, indicate that the contribution of nonther- 
mal processes (UV pumping) to the II2 excitation is rather 
moderate. 
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The excitation conditions in NGC 1275 can be investi- 
gated in more detail by examining the ratio of two ortho 
or para lines and put them in context. The remaining un- 
certainty introduced by the intrinsic ortho-para ratio of 



the gas clouds can thus be avoided (Mouri 1994). Locat- 
ing these ratios in a ortho-para line ratio plane provides 
an independent guide for distinguishing between the rela- 
tive contribution of X-ray heating, shocks, and UV pump- 
ing excitation of H2 gas in starburst and AGN galaxies 
(Fig. As non thermal processes (e.g. UV pumping in 
low density gas) excite high vibrational states more easily 
than the thermal process, the ratio of two ortho transitions 
such as Ro = H2 2-lS(l)/H2 1-0 S(l) gives an ortho-para 
balance independent measure of the relative contribution 
of thermal and non-thermal processes. The same holds for 
Rp = H2 1-0S(2)/H2 1-0 S(0). The expected value of Rq 
in the case of pure non thermal excitation under nearly all 
realistic conditions is Rq = 0.56 (Black & van Dishoeck 



1987), while thermal UV excitation gives a value of Rq 

|1987D . From Table | we find 
14.4/10.0 = 1.4, which 



al. 



near (Sternberg et 
Ro = 7.7/41.1 ^ 0.19 and Rp 



again, as we noted earlier, puts NGC 1275 in the Ro-Rp 
plane close to the location of NGC 6240 (Fig. ||). Starburst 
galaxies like NGC 253 and galactic photodissociation re- 
gions generally show Rp < 1.0. Following this method one 
can also define a ratio, Rq — H2 1-0S(3)/H2 1-0 S(l), of 
two ortho lines which is again independent of the initial 
ortho-para ratio in the molecular clouds. When placed in 
this diagram (Fig. NGC 1275 falls near other powerful 
AGNs like NGC 1068. 

The interpretation of these diagrams tells us that the 
excitation is mostly "thermal" , with a small component 
appearing as "nonthermal emission" . The gas must then 
be quite dense, with n^p >, lO^cm"'^. The placement of 
NGC 1275 near other powerful AGN and merger galaxies 
suggests that both X-ray and shock mechanisms may be 
at work creating the powerful H2 emission. Hence, based 
on the spatial extent of the H2 emission and its tempera- 
ture, we propose that there are two distinct sources of H2 
excitation: the AGN nucleus, which excites the vast ma- 
jority of the radiation, and which produces the compact 
emission core seen in our maps, and the extended region 
of moderate star formation which excites the highly ex- 
tended, weak H2 envelope. 



J^.'S. Hi 2-lS(3) line strength 

The column density of the H2 2-1S(3) line is lower by 
a factor of 3 compared to what could be expected from 
weighted extrapolation of the 2-1 S(l) and 2-1 S(2) tran- 
sitions. Even if the ortho/para ratio within the 2-1 transi- 
tions is not close to 3, which it was in the 1-0 transitions, 
N, 



for blending the line with a strong absorption can be ex- 
cluded for this part of the spectrum and the atmosphere, 
we argue that this weakness of the H2 2-1 S(3) line is real. 
Inspecting Fig. ^reveals that the 2-1 S(l) line is also weak 
outside the nucleus at a radius of « 1 kpc. It would be vis- 
ible if it was 2.5 time stronger or about half the strength 
of the 1-0 S(2) line, which is clearly present. 

The only other galaxy, were a weak H2 2-1 S(3) line has 
been reported so far, is NGC 6240 (Lester et al. 1988), al- 
though only as an upper limit. From their Table 3 one 
can estimate that the column density of 2-1 S(3) line in 
NGC 6240 is at least 3.5 times weaker than expected from 
extrapolating other data, which is in good agreement with 
our finding and suggest that there may be a common 
mechanism at work. The column densities in NGC 6240 
are only a factor of 2 higher compared with NGC 1275 an 
their derived temperature range is almost identical with 
Fig. ||. It should be noted, however, that the excitation 
conditio ns in both galaxies are not identical. Draine & 
Woods ( 1990| ) have shown that in addition to X-ray heat- 
ing, UV pumping plays a major role in NGC 6240. This 
is already indic ated in Fig. || as well as and in Fig. 4 
of Lester et al. ( 1988| ): The large offset between the line 



through the 1-0 transitions and the 2-1 data is diffe rent 
from NGC 1275 (Fig. |). Recently, Sugai et al. ( |1997D re- 
ported that they were unable to confirm the weakness of 
the H2 2-1 S(3) line in NGC 6240 with R = 580 grism spec- 
troscopy. Better quality data at higher spectral resolution 
is certainly required to clarify the situation. 

The most likely mechanism for the weakness of the 
Il2 2-lS(3) line is resonant fluorescent excitation of the 
upper le vel {u = 2, J= 5) by Lya photons (Black & van 
Dishoeck 1987). The Lya photons needed for this pumping 
process can easily be provided in an environment where 
a notable fraction of the II2 is excited by X-ray heating 
(Black & van Dishoeck 1987). Although detailed models 
of resonant fluorescent excitation of the = 2, J = 5 level 
does not exist, the weakness of the II2 2-1 S(3) line by 
itself already demonstrates that at least in the nuclear re- 
gion the excitation of the molecular hydrogen in NGC 1275 
is dominated by the central AGN. This is in agreement 
with our previous finding that the excitation environment 
in the nuclear region of NGC 1275 is much more dom- 
inated by X-ray heating than in NGC 6240 (Draine & 



Woods 199C). NGC 1275 may be the first galaxy for which 



a more complete model on resonant fiuorescent excitation 
can be computed. Finding other galaxies with such a weak 
H2 2-1 S(3) line (e.g. in NGC 1068) would be important 
for constraining the process, which may provide an ad- 
ditional tool for understanding the radiative interaction 
between AGN and circumnuclear molecular material. 



„_i(2-l S(3)) is still a factor of > 2.5 lower compared to _ . . 

)n io/i(( J • nnnn tri: T ^- Br7emiSS10n 



NpQi(2-l S(l)) and using Tex = 2600 K for the two ortho 
states. Inspecting the original data and the spectrum in 
Fig. |l| did not show any abnormalities. Since the chance 



The Br7 emission from NGC 1275 is unusual in that it 
is highly asymmetric both in its line profile (Fig. |l|), and 
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Fig. 8. Upper level column 
density per mode of observed 
H2 lines against upper level 
energy in K. A slope in this 
diagram is equivalent to an 
excitation temperature. 



in its spatial distribution (Fig. ^0|). The Br7 emission is 
extended along roughly the same line (PA =160°) as the 
powerful radio jets (Pedlar et al. 199C), which suggests 





Fig. 9. Excitation diagrams showing the relative contri- 
butions of thermal and nonthermal excitation sources to 
H2 emission (styled after Mouri 1994). The position of 
a source in these diagrams is independent of the intrin- 
sic ortho/para ratio in the molecular clouds. In both dia- 
grams the regions occupied by shocks (Brand et al. 1989), 
nonthermal (Black & van Dishoeck 1987), thermal UV 
(Ster nberg & Dalgarno 1989 ), and X-ray (Lepp & McCray 



1985, Draine & Woods 1990) excitation are outlined. Data 



of NGC 6240 and NGC 1068 are taken from Lester et al. 



(1988) and Oliva & Moorwood ( |1990D . The location of 
NGC 1275 in these diagrams speaks for a largely thermal 
excitation of > 2000 K. 



that the mechanism of its ionization is more strongly in- 
fluenced by the central engine than by the surrounding 
stellar population. This result correlates well with Ha and 
[Nil] maps with visib le band imaging spectroscopy (Fer- 
ruit & Pecontal 1994). In particular, the position of the 
our maps correlates with an Ha and [Nil] 
|ak located « T.'2 N, 0'.'7 W of the nucleus, 
blished Ha maps, we estimate the Ha flux in 
to be 1.3 X 10"^'' W m~^ arcsec"^. Based on 
the observed [SII] line ratios and assuming a temperature 
of Te = 10** K, (Ferruit & Pecontal |1994| ) derive a mean 
electron density of ne — 40 for this region, and a typical 
linewidth of 340 km s^^. 

The intrinsic Br7 linewidth at the nucleus of 980 ± 30 
km s~^ is unusually large for a pure star forming source. 
However, the width is not unexpected for emission from 
NLRs associated with AGNs and we therfore safely as- 
sume here that the NLR contributes a substatial fraction 
to the Br7 emission. The larger linewidth of Br7 compared 
to Ha is probably a differential extinction effect. The emis- 
sion was too faint to derive accurate Br7 line profiles in re- 
gions 1" to the NW and SE of the nucleus in order to test 
whether the excess width and shape of the Br7 line is due 
to expansion. The total Br7 emission within the central 
1 kpc is 9.7x 10~^^ ergs""'^ cm~^, which is a total luminos- 
ity of 1.5 X 10^ Lq. This corresponds to an ionization rate 
of 3.2 X 10^3 g-i or 2.0 x 10^ Lq for Tg = 7500 K, and can 
be compared with the value of 1.1 x 10^* s~^ inferred from 
M82 ove r a si milar region from H53a recombination lines 
(Puxley |l99l| ) . From the ionization rate one can infer the 
total Lyman continuum luminosity ^Ijjq if '^^^ assumes 
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burst/Seyfert indicators than it is in most active galax- 
ies. For example, in the central Ikpc of NGC 1275, we 
find [Fell](A1.644)/Br7 = 6.3, whereas in roughly simi- 
lar sized regions in the merger driven starburst NGC 3256 
[Feli]/Br7 = 1.55 and in the starburst/Seyfert systems 
NGC 4945 [Feii]/Br7 = 1.4 (Moorwood & Oliva ^99§ 
and NGC 1068 [Fei i]/Br 7 = 2.5 (Thatte et al., in prepa- 
ration; Blietz et al. 1994). While Br7 originates predomi- 




Fig. 10. Line map of Br7 showing clear asymmetrical 
extension. Solid line marks the radio jet position angle 
of 160°. Contours are at 10, 20, 30. . . % of peak flux of 
5.3 X 10~^^ erg s^-^ cm^^ arcsec"^. 



OB type stars as its source. LLy^ — 2.0 x 10^ Lq yields 



nantly from photoionization of H ii regions by OB stars, 
[Fell] has an additional source. Most of the gas phase 
Fe in the local ISM condenses onto grains, so that the 
free Fe is often depleted by factors of 100 - 1000. If we 
presume that the ionization characteristics of the central 
source in NGC 1275 are similar to those of NGC 3256, 
NGC 4945, and NGC 1068, then the u nusua l [Feii]/Br7 
and [Feii]/Pa/3 ratios (e.g. Rudy et al. 1993 ) would sug- 
gest that the iron depletion is ^ 10 in NGC 1275, less 
severe than in most starburst and active galaxies. This 
could be due to an ambient radiation field which heats 
the dust grains and tends to evaporate the Fe off of them. 
Such a scenario is supported by the relatively high de- 
rived dust temperature within the central 300 pc (e.g., 
Sect. ^^). Alternatively, the fact that the ionization state 
of the gas must allow for plentiful Fe+, which can then 
be further ionized by photons with E > 16.2 eV, requires 
that the ionizing source of radiation must be relatively 
cool, in order to produce low ionization radiation which 
would produce the high [Feii]/Br7 line ratios. However, 
fast J-shocks and X-ray heating can also produce copious 



about 1.5% of the Lpjj^ of NGC 1275 (1.5 x lO" L0) or [Fell] emission (Blietz et al. |1994D. In the latter case, the 



^FIr/ ^Lyc ^ ''5- ^Lyc 



2 X 10^ Lq is a lower limit to 



the Lyc-luminosity in the IRAS beam. The total bolomet- 
ric luminosity of such an evolving stellar cluster is typically 
about 10 - 20 times higher than the Lyman continuum lu- 
minosity for typical assumptions on the properties of such 
a region if it were an active star forming region powered 
by massive stars (age ~ a few 10'' y, Mj^^^j. = 0. 1 Mq , 



Mupper = 120 M0, see Fig. 10 in Krabbe et al. |1994D 



With that ratio Lg^j^^^^ ~ 3.8 x IO^'^Lq, which is about 
25% of Lpjj^. Taken as face values these numbers are not 
very consistent with star formation but much more with 
what one gets empirically, using above formulas for AGNs 



(Genzel et al. 1998). 



6. The iron emission 

6.1. Strength and morphology of the [Fell] emission 

The [Fell] 1.644 fj,m emission from NGC 1275 is very 
strong (Fig. |^) and centrally concentrated (Fig. |ll|). Other 
prominent H and K band iron lines include the [Fe 11] 
1.534 /Ltm, [Feiii] 2.2178 /xm, and [Feiii] 2.2420 fim, 
in a near blend with Il2 2-lS(l). The flux of the [Fell] 
1.644 /im (Table ^) is nearly as large as the [Fell] 
1.257 /zm in an 8" aperture by (Rudy et al. |l993D . The 
[Fell] emission is much stronger compared to other star- 



X-rays penetrate deep into molecular clouds where they 
create large reservoirs of partially ionized gas. 

6.1.1. Profile and excitation of the [Fell] emission 

The profile of the [Fell] 1.644 /im line cannot be well 
modeled as a single Gaussian. This is typical of Seyfert 
galaxies, which often show narrow and broad components 
of the same emission feature. We therefore fit the observed 
profile using two Gaussians, a broad component of width 
1340 ± 120 km s~^ and relative amplitude of 0.26, and a 
narrow component of of width 280 ± 20 km s~^. This two 
component fit (Fig. |l^, is excellent, and suggests that 
w 50% of the [Fell] emission arises in the same region as 
the Br7, showing a similar linewidth, while the other half 
of the [Fe 11] arises in the same region as the II2, where the 
linewidths are again similar. Hence we propose — as in the 
case of the H2 lines (Sect. 4.2) — that the active nucleus of 
NGC 1275 is responsible for much of the excitation leading 
to [Fe 11] emission. Comparing this result with Blietz et al. 
(1994), who have found an excellent correlation between 
the [Fell] emission and the NLR in NGC 1068 suggests, 
that we might have a similar case here but unresolved 
and at much greater distance. 

The detailed line ratios expected from excitation via a 
central AGN have been modeled by HoUenbach & Maloney 
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Fig. 11. Line map of the Fell 1.644 /im line showing faint 
extensions beyond the immediate nuclear area. Contours 
are at 1.25, 5, 10, 20, 30. . . % of peak flux of 3.0 x lO^^^ 



erg s ^ cm ^ arcsec ^. 



( |1996D . In their model, the X-ray energy absorbed per unit 
time per hydrogen nucleus in a cloud at distance lOOrioo 
pc is Hx 7 X 10-22 L44rj"ooN^2^ ergs-\ where 10'"'L44 
erg is the AGN luminosity and IO22N22 cm-^ is the 
total column density. The heating rate is w 0.3-0.4 Hj^, 
and the cooling is due to line emission from many species, 
but of particular diagnostic interest are the H2 1-0 S(l) 
and [Fell] 1.64 /im lines. Using our value of N j^^^ = 

H2 

3 X 10^^ cm~^, a molecular abundance of H/H2 ^ 10^, and 



?a 2 X 10 erg s (Prieto 1996 ), and integrating the 



HoUenbach & Maloney results over our resolution beam 
of 300 pc, we can account for the observed [Fell] to H2 1- 
0S(1) line ratio. For all radii larger than our beam, the 
predicted ratio is « 1 : 1. This is close to our measured 
value of 1.5:1 in a 3"aperture (Table ||) since the predicted 
ratio depends on the presumed depletion of Fe, which we 
already know to be < 10, and which the model presumes 
to be « 30 (Hollenbach, private communication). We thus 
conclude that the [Fell] emission as well as [Feli]/H2 1- 
0S(1) can be explained by an X-radiation field typically 
emitted by AGN. 



7. The COsp index 

The COsp index of the underlying stellar population can 
be diluted by continuum emission not originating in stars. 
If the non-stellar continuum is x times stronger than the 



Fig. 12. Profile of the Fell 1.644 /xm line showing mul- 
ticomponent fit. The solid dark line is the data, to which 
the best single component fit has a full width of 544 km 
s~^. The two lower intensity profiles are the broad and 
narrow component required to fit the profile well. The 
dotted line which closely follows the data is the sum of 
the broad (1345 ± 120 km s'l) and narrow (280 ± 20 km 
s"^) components. 



stellar continuum, the COsp index will be diminished ac- 
cording to 



-2.5 log 

-0.4COsp(0) 



Sq+x 
1 + X 



(2) 



where So = 10"" '""^spi^^J ^j^g 

average of the undiluted 
rectified spectrum between 2.31 and 2.40 /im and COsp 
is the spectroscopic index. Details are given in Doyon 
et al. (1994). According to our model fitting, the non- 
stellar continuum in the central 1" contributes ^ 90% 
of the emission, reducing the observed COsp from its in- 
trinsic value of w 0.3 to w 0.03, and and giving rise to 
the COsp hole at the nucleus of NGC 1275 (Fig. ^). The 
range for the allowed values of the intrinsic COsp is in- 
dicated in Fig. The nuclear stellar cluster may have a 
higher CO index than initially assumed, implying that the 
K supergiants with COsp ~ 0.3 instead of giants may be 



the dominating population (Doyon et al. 1994). Including 
this finding into a reiteration on the continuum decom- 
position will not have any effect on the results, since the 
temperature of these stellar types are not very different. 

Outside of the hole, the COsp of 0.17 is typical of 
E/SO galaxies (Frogel et al. 1978), confirming that the 



stellar population of NGC 1275 is quite normal. Further- 
more COgp = 0.17 corresponds to a K5III star (Doyon 
et al. 1994 appendix A), which corroborates our spectral 
synthesis results for regions at R <L 3". 
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Given that the hght is dominated by K giants or super- 
giants, we try to estimate the total nuclear stellar mass. 
Using our result that the emission enclosed in the cen- 
tral 3" (Ikpc) of NGC 1275 is > 90% due to non-stellar 
sources, we correct the observed K magnitude of the stel- 
lar component only of mj^ = 11.8 [in agreement with 
other recent values mj^ — 11.9 (Forbes et al. 1992)] to 
mj^ = 14.3. Converting to absolute visual magnitudes 

yields M^SHI (jv^qq 1275) = (V-K) + mj^ - 25-Aj^ - 5 
log T>Mpc = "16.8 using V-K= 3.24, Aj^ - 0.2, and 
D = 68.8 Mpc. The mass of a K giant is 1.2 Mq, and its 



M 



K5III 



y - -0.2 (Lang |l992| ), which implies 4.6 x 10^ Mq 
of K5III stars in the central 3". 

Assuming a Salpeter IMF with P(M) oc M-i-^'^ ^^y[^ 
the mass fraction of the stars in the range 1.0 to 1.5 
Mq (roughly the range of stellar masses which will be- 
come K5 III stars), for a population bounded by 0.1 Mq 



-0.35 



1.5 



-0.35 



)/(0.1 



-0.35 



and 100 Mq stars is (1.0 
200-0.35) = 0.065, so that K5III stars make up about 
7% of the mass of a star cluster. Hence the total mass of 
stars in the central kpc is thus 7.1 x 10^ Mq. If the pop- 
ulation is dominated by K supergiants, the total stellar 
mass comes out to be the same value (mass range = 12 - 
14 Mq M^^I = -5.0). 

We estimate the core radius ro of the nuclear stel- 
lar cluster to be smaller than 500 pc. Using this value 
as an upper limit and knowing the stellar mass of the 
cluster, we can use po = 9(TQ/47rGrg (Eckart et al. 



1993) to estimate the expected velocity dispersion: a — 
V -^cluster ^/^'''o ~ 8.4 km/s. Even if the CO absorption 
had been detected on the nucleus, such a small dispersion 
had been beyond the capabilities of our instrument. 

We conclude that the absence of CO bandheads in 
the nuclear emission indicates that those bandheads are 
heavily diluted by the emission of hot dust and power-law 
emission. 



8. Discussion 

Our spectral synthesis (Sect.|3[) shows that inside of the 
central 300 pc, the NIR continuum emission is dominated 
by the AGN and hot dust components. Together these are 
10 times more powerful than the stellar emission from 
the same region. This result is confirmed by the observed 
factor of ^ 10 continuum COsp dilution (Sect.Q) at the 
nucleus. The nature of the nuclear power source likely is a 
massive black hole (MBH) with a surrounding molecular 
gas cloud torus or disk. The lack of strong CO absorption 
bands indicates that massive star formation probably does 
not play a very important role in the nucleus. 

The intensity and morphology of the line emission sug- 
gest that more than one mechanism is at work exciting 
the line emission. In addition to the AGN, shocks are also 
likely to play a large role. Since the J band [Fe ii] lines are 
unusually strong, even for an AGN, it requires a substan- 




Fig. 13. The azimuthally averaged COsp index in 
NGC 1275. The apparent decrease in COsp toward the 
center is an artifact of continuum dilution due to hot dust 
and the scattered light of the AGN. 



tial boost of the II2 emission to keep the [Feii]/Il2 1-0 S(l) 
ratio near unity. This is indeed observed in NGC 1275, 
which has unusually strong H2 emission compared to other 
Seyferts. Furthermore, the spatial extension of the H2 



emission (Sect. 4.1) suggests that it may be driven by 



shocks from a recent or ongoing merger event. Such a pic- 
ture is supported by not only the o bserved isophotal shell 
structure (Hernquist & Weil 1992 ), but als o by the ex- 
tinction analysis of Norgaard-Nielsen et al. ( 1993| ), which 
shows that the HV system is at least partially moving 
through the main galaxy rather than completely in front 
of it. Hence the interaction of the HV and LV systems 
could give rise to the required shocks which boost the 
H2 emission. This process is postulated to occur in the 
merging galaxy NGC 6240 (van der Werf et al. |1993D. 



Among AGNs, strong [Fell] emission is not unique to 
NGC 1275. In fact there are a large number of AGNs with 
elevated [Feii] emission (among them many IRAS lumi- 
nous and starburst galaxies such as MRK231, MRK507, 
and PHL1092), which can be attributed to violent star 
formation in a metal rich environment (Fillipenko & Ter- 
levich 1992| , Lipari et al. 1993). In the merger induced H2 
boosting scenario, the unusually high level of [Fell] emis- 
sion observed from NGC 1275 would likely be due to a 
different effect, namely the dissociation of Fe from inter- 
stellar grains via shock h eating . This could produce an 
iron rich gas phase (Sect. |6TT| ), which would then give 
rise to the very strong observed [Fe 11] emission. 
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because 200 Mq of star formation would be 20 times more 
powerful than the prototypical starburst galaxy, M 82, and 



Fig. 14. The black lines in the graph represent solutions of 
equation (H). The inner area of the triangle illustrate the 
range of allowed values of the measured CO spectroscopic 
index COsp(x) in NGC 1275. Constrains are the limits set 
by the error of the index itself, the range of likely values 
for the spectroscopic index COsp(O) in a nondilute d ste llar 
cluster, and the error in x (see Fig. |l^ and Sect. 3.2,) 



The fact that the Br7 emission is generally extended 
suggests the presence of young stars and thus extended 
recent star formation. If the slight extension of the con- 
tinuum is significant, it suggests increased star formation 
along the jets, perhaps due to shocks, caused by interac- 
tions with the surrounding medium. Part of the Br7 emis- 
sion along the radio jet may be influenced by the radiation 
from the nucleus and thus not created by star formation. 
The extended component of hot molecular hydrogen is 
probably shock excited as a result of the star formation. 

The existence (or non-existence) of "cooling flows" — 
an inflow of cool gas into a galaxy from its halo of hot 
x-ray emitting gas — is currently a topic of intense de- 
bate. The evidence for them is that some galaxies which 
are surrounded by a halo of hot. X-ray emitting gas show 
a dramatic central drop in the gas temperature. As the 
primary cooling mechanism of this gas is bremsstrahlung 
(emission cx n^), and as the central temperature drop is 
accompanied by a brightness increase, the gas must be 
both cooler and denser (see, e.g. Sarazin et al. 1992 for a 
physical introduction and review) . If the gas is cooler and 
denser, then it must fall toward the center of the galaxy. 
In the case of NGC 1275, the predicted infall amounts to 
« 200 M0yr"i (Fabian et al. |1981D , or 2 x lO^^ 
in 10^° years. The evidence against cooling flows is that 
if such enormous quantities of gas really are falling in, 
they must be detectable in some way, yet careful searches 
have not revealed the gas in any of its expected states. We 
know that the gas cannot form stars with a normal IMF, 



easily detected via color variations (Cardiel et al. 1995). 

Our II2 data provide several evidences against the cool- 
ing flow hypothosis. 

i) The spatial distribution of the hot molecular hydro- 
gen strongly peaks at the nucleus of NGC 1275 on a scale 
<300pc. If the H2 emission was due to the cooling flow, 
one would have expexted a much more extended distribu- 
tion since the gas has to fall through all of the galaxy. The 
weakness of the H2 2-1 S(3) line (see Sect. 4.3) emphasises 
the close relation of the hot H2 with the nucleus. 

ii) The temperature variation within the hot H2 is 
small. Between K and about 20000 K upper level energy 
the temperature only varies by about 1500 K. A cooling 
hot gas cloud is likely to show a much larger span of tem- 
perature compared to what our data show. 

in) The temperature of the hot H2 is rather low: be- 
tween 1500 K and 3000 K. There is no evidence of a very 
hot ( > 3000 K) II2 component which can be expected 
for a cooling hot gas mass. 

iv The line emission of the hot H2 can be explained 
consistently with mostly thermal excitation (e.g. in a post- 
shock environment) with a smaller component of X-ray 
heating. A "cooling flow" scenario does not need to be in- 
voked to explain the line emission of molecular hydrogen 
in NGC 1275. 

Several groups have argued that the gas cannot have 
collected in the form of molecular hydrogen (H2), because 
searches for CO emission associated with H2 have revealed 
only very limited quantities of H2 gas, although NGC 1275 
is unusual in that it shows more CO than most cooling flow 
galaxies (Lazareff et al. 1989| , Inoue et al 1996, O'Dea et 
al. 1994). However, since the temperature regime is very 
different from the cold molecular clouds, the CO H2 
conversion factor is certainly different and may even not 
be applicable in this context. 

The gas could have been collected as atomic hydrogen. 
However, only very limited quantities have been detected, 
certainly much less than would be expected from an in- 
fall of several hundred solar masses per year over cosmic 
timescales (Haschick et al. 19821 ). 

One of the last places in which the purported inflow 
could be hiding is in stars with a very unusual kind of 
IMF, one that is biased toward producing low mass stars, 
so that the star formation rate as determined by the total 
number of ionizing photons (these are produced only by 
massive stars) does not reflect the total amount of mass in- 
volved in the star formation process. A population of such 
low mass stars should have a somewhat lower COsp index 
than a normal population, because it will have fewer lumi- 
nous giants and supergiants which are the largest contrib- 
utors to the CO absorption lines. Our NIR spectroscopy 
allows us to analyze the stellar population of regions 1 kpc 
from the nucleus via the COsp index, where we flnd a 
COsp = 0.17. This can be compared with theoretical pre- 
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dictions for a population of very low mass stars (see Fig. 
20 in Kroupa & Gilmore 1994), whose expected value of 
COsp is between 0.07 and 0.16 for cooling flow rates of 
158 Mq yr~^. For higher mass deposition rates, the ex- 
pected COsp decreases, so that for a rate of 316 Mq yr~^ 
the expected COsp lies between 0.0 and 0.1. Hence we 
cannot completely eliminate the existence of a low mass 
star population for a cooling flow of 200 Mq yr~^ as is 1985) 
postulated for NGC 1275, but the observed COsp index is 
at the extreme limit of theoretically acceptable values. 



9. Summary 

We observed the central 1 kpc of NGC 1275 with the 
MPE 3D imaging spectrograph. We find that the cen- 
tral w 300 pc are strongly dominated by emission from 
an AGN and from hot dust: There is no evidence of a 
nuclear stellar continuum. In the regions at a distance of 
~ 1 kpc from the nucleus, the situation is completely re- 
versed, and the emission is totally dominated by a stellar 
population whose spectra is consistent with a normal, old 
(> 10^ — 10^ y) population. In this region, we find no 
evidence to support the thesis that substantial quantities 
of infalling "cooling flow" gas has been transformed into 
stars with an IMF biased toward low masses, although we 
cannot strictly rule out a low mass population because the 
purported flow may be disturbed in the central 10 kpc of 
the galaxy. Within the central ~ 500 pc, the emission line 
ratios of many important diagnostic lines from the NIR 
regime are consistent with excitation from a combination 
of AGN emission and shock excitation. The shocks may 
be due to previous mergers, or to the early phases of a 
merger between the HV and LV systems. 

Acknowledgements. We thank the 3D team for their enthusi- 
astic support during the observations. 



to obtain the synthetic spectrum. The extinction has two 
parts, one due to dust within our Galaxy [grx), and one 
due to dust within the source galaxy {t\). In our model, 
the latter can have either a simple screen geometry or a 
mixed stars/gas geometry. In all cases, the extinction is 
taken to vary as A~^'^^ (Landini et al. 1984), and is fixed 
to be Apj = 0.112 Ay at A 2.2 ^m (Rieke & Lebofsky 



When comparing l'^^ with the observed spectrum, 



rOBS 



we avoid emission lines due to non-continuum 



sources (e.g. shocked ISM gas, photoionized gas clouds), 
and define our goodness of fit parameter, x^, such that 



1 



(,SYN 



0BS)2 



N -I 



(A4) 



X-line 



where is the number of measured spectral points 
which contain no line emission or absorption from non- 
continuum sources and a\ is the local noise. 

To model the stellar component, we have resampled 
the R« 3000 normalized K band spectra of Kleinmann 



HaU ( |1986| ) at the 3D resolution of R = 6 70 ont o our 

p^96| ). 



0.001 /im pix~ data format (Forster et al. 1996|), and 
then multiplied them by a Rayleigh-Jeans power law of 
^■^'^ (Doyon et al. 1994) corresponding to an AOV star 



A 

to reconstruct the true stellar spectrum. 

The intrinsic spectrum of the hot dust of temperature 
Tj, is taken to be a black body multiplied by a power- 

law emissivity of the form Cfj^st A~ . Hence, = 
BxiT^X'^ , where Ba(T^) is the black body function. 

The behavior of free-free emission is given by /j, oc 
g-hiz/kT -^yj^gj-g ^ jg ^]2e plasma temperature (normally 
> 5000 K). We convert to wavelength units, and find Ix oc 
A^^e"^'^/^'^^ which we adopt as our functional form for 
free-free emission. 



Appendix A: Spectral Fitting of the continuum 
emission 

The total intrinsic continuum emission at each wavelength 
A is 

^l^^^E'^i^i ^^^'^ E'^i-i- (Al) 

i i 

In active and starburst galaxies, we identify four sources 
of continuum emission: stellar (J^'l'A^), AGN power- 

law(/P) and hot dust [I^^], and free-free [if^]. The 
intrinsic emission of 

rINT_^ rSTAR ,^ jP , ^. rHD , ^ rFF/A9^ 
is then extincted according to 



.SYN ^ f/PTg-gr^g-r^^ Screen /^o^ 

^ U{^^e-9-A(l-e---)/TA, Mix. ^ ' 
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